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ABSTRACT 



This thesis develops an aggregate demand inventory model for 
repairable items. It uses minimization of wholesale stock 
investment level as the objective function subject to a 
given mean supply response time (MSRT) goal. Also addressed 
are annual budget constraints encountered by Inventory 
Control Points (ICPs) as well as a constraint placed on the 
total wholesale investment level which is implied by ceil- 
ings on the Navy Stock Fund (NSF). Preliminary parametric 
analyses of the model showed that the wholesale stock 
investment levels increase at a decreasing rate as repair 
induction batch sizes are increased and attainable MSRT 
values decrease exponentially as investment levels increase. 
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I. INTRODUCTION 



A . BACKGROUND 

The Department of Defense, and specifically the Navy, is 
facing one of its greatest challenges. As increasing defense 
dollars are authorized and appropriated for both old and new 
weapon systems and their support, many in Congress and the 
general public are demanding that maximum measurable benefit 
be received from each dollar expended/invested . In partic- 
ular, the military services are being asked to demonstrate 
an improvement in readiness and sustainability consistent 
with the increases in investment. These increases in invest- 
ment have come in many resource areas ranging from manpower 
to weapon systems. 

One such resource area is the investment in spare parts 
which is required to ensure the operational readiness of the 
Navy. As increasingly complex weapon systems are introduced 
into the Navy, the availability of spare parts to support 
these systems becomes all the more critical. If complex and 
expensive systems are not operational when they are needed, 
then all funds invested in their research, development, 
production and introduction into the Navy will have been 
basically for naught. Therefore, the proper management of 
spare parts is of prime importance to the Navy in accom- 
plishing its mission. 
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The rapid expansion the Navy has experienced in recent 
years combined with significant advances in technology have 
resulted in much more sophisticated equipment. This increase 
in sophistication has increased the complexity of the Navy 
Supply System's efforts to maintain sufficient stocks of 
replacement components and repair parts. In particular, more 
and more attention is being focused on the management of 
depot-level repairable spare parts. Items designated as 
depot-level repairables or DLRs are those items which must 
be removed from their weapon system and returned to a desig- 
nated overhaul point (DOP) for repair when they fail. 

The uniqueness of the military requires that its objec- 
tive in inventory management differ from that in the private 
sector. In the private sector the objective utilized in 
inventory models is cost minimization due to the profit 
motive of private firms. However, in the Navy, the Chief of 
Naval Operations (CNO) has directed that Supply Material 
Availability (SMA) be the measure of effectiveness at the 
wholesale level [Ref. 1]. The wholesale level contains 
back-up inventories which may be requisitioned by any 
customer worldwide. 

Although SMA has been specified as the measure of effec- 
tiveness, the Navy Inventory Control Points (ICPs) which are 
assigned responsibility for wholesale management of 
repairable items continue to use inventory models based 
primarily on cost minimization and relate SMA to backorder 
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cost. This is a consequence of the DOD Instruction 4140.39, 
"Procurement Cycles and Safety Levels of Supply for 
Secondary Items" which specifies a consumable model of mini- 
mization of annual ordering and holding costs subject to a 
constraint on time-weighted, essentiality-weighted units 
short. No instruction has been written for repairables but 
the Navy assumes it would be of a comparable form. However, 
alternative objective functions can be considered which 
incorporate SMA or time-weighted units short. The recent 
wholesale initial provisioning model developed by Richards 
and McMasters [Ref. 2] and approved for Navy implementation 
in December 1984 had as its objective the minimization of 
mean supply response time (MSRT). MSRT is directly related 
to time-weighted units short. As a consequence, MSRT was 
also chosen by Apple as the objective function for replen- 
ishment of repairable items at the wholesale level [Ref. 3]. 
His model serves as a foundation for this thesis. 

B. OBJECTIVES 

Apple [Ref. 3] proposed an inventory model for manage- 
ment of repairable items at the wholesale level which is 
readiness - vice cost - oriented. However, the Navy cannot 
focus on readiness alone, because it must operate within 
various budgets which are authorized by Congress. These 
budgetary constraints, both annual procurement/ repair 
budgets and long-term budgets implied by Congressional ceil- 
ings on the Navy Stock Fund (NSF), must be considered in any 
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model designed for inventory management in the military. 
Therefore, the first objective of this thesis is to expand 
upon Apple's model by converting it to an aggregate- demand 
model for repairables management which seeks to minimize 
wholesale stock investment subject to budgetary constraints 
and a mean supply response time goal. The second objective 
is to conduct parametric analyses of the model to ascertain 
how it will perform under a variety of parametric changes. 

C . PREVIEW 

A brief discussion of the repairables system is provided 
in Chapter II. It is followed by a detailed review of the 
multi- echelon wholesale model proposed by Apple which has 
minimization of mean supply response time (MSRT) as its 
objective. Chapter III presents the conversion of the Apple 
multi- echelon model to an aggregate-demand model which can 
be used at the ICP level with the current UICP demand fore- 
casting program. It also proposes minimization of wholesale 
stock investment levels as the objective. Finally, it 
formulates the annual budget constraints on procurement and 
repair which are encountered by the ICPs . The model devel- 
oped in Chapter III is analyzed parametrically in Chapter IV 
and the results and conclusions of the analyses are 
discussed. Chapter V provides a brief summary and some 
final observations and recommendations of areas for further 
study of the aggregate-demand model. 
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II. MEAN SUPPLY RESPONSE TIME REPAIRABLES MODEL 



A. THE REPAIRABLES SYSTEM 

An item of supply is designated as a repairable if it 
can be repaired faster and/or less expensively than it can 
be procured. The management of repairables begins with the 
designation of an item as a repairable during Weapons System 
Acquisition and continues with the procurement process and 
the repair cycle. This complete management system is called 
the repairables system. 

Repair is accomplished at one of three maintenance 
levels: (1) the organizational or lowest level (i.e., a 
ship); (2) the intermediate level (i.e., a tender, carrier, 
or a shore Intermediate Maintenance Activity); (3) the depot 
level (i.e.. Naval Shipyard, Industrial Naval Air Rework 
Facility or a commercial repair activity). Figure 2.1 
depicts the repair cycle for a depot level repairable (DLR) . 

The process begins when the ship or customer registers a 
demand for a DLR at the nearest stock point (NSC) (which is 
the point of entry) . If the item is available directly from 
that NSC, it is issued to the customer. The demand for an 
item that is not available is referred to the inventory 
manager at the ICP (SPCC) who must either refer the requisi- 
tion to another stock point holding the item or record the 
requisition as a backorder against stock that is due-in. 
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Figure 2.1 Repairables Cycle. 

Once the demand for the repairable is satisfied, the 
inventory manager is concerned with the customer getting the 
not - ready- for- issue (NRFI) carcass into the repair cycle. 
The repair process for a DLR begins when the NRFI is shipped 
to a specified NSC where it is held until a predetermined 
quantity, R, (based upon existing inventory management 
policy) is available for induction to the Designated 
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Overhaul Point (DOP) for repair. Once repairs are complete, 
the DOP ships the ready-f or-issue (RFI) unit to the NSC 
designated by the inventory manager. 

Since there is some attrition of the repairables in this 
closed loop system due to items being lost or being beyond 
economic repair, the inventory manager must ensure that a 
fixed level of total units, SW (both RFI and NRFI), are 
available within the wholesale system. This is accomplished 
by procurring a quantity, Q, from a manufacturer whenever 
the wholesale system stock falls below a predetermined 
reorder point, ROP = SW - Q. The manufacturer ships these 
units to the NSCs as directed by the inventory manager. 

Let the items defined below represent the average times 
it takes for the described events to occur: 

Tl: carcass turn-in time; i. e. the time it takes for 
a carcass to be received at the collection point 
(Naval Supply Center - NSC) after a demand has 
been registered (this includes shipboard turn-in 
time and shipping time); 

T2: shipping time for a carcass from the NSC to the 
DOP; 

T3: shipping time for an RFI unit from the DOP or a 
manufacturer to the NSC; 

T4: shipping time for an RFI unit from the NSC to a 
ship ; 

T5: time required for the ICP to determine that a 
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carcass will not be returned to the system; 

RTAT : time required for the DOP to repair an item or a 
batch and return the batch to RFI condition. 

ALT: administrative lead time required by the ICP to 
prepare a purchase order or contract and the 
ordering data to purchase a replacement item; 

PLT : production lead time required by the manufacturer 
to manufacture the quantity of an item being 
purchased . 

B. APPLE'S THESIS 

Apple [Ref. 3] provides a detailed overview of the 

repairables system and its importance to the Department of 
Defense readiness, the roles of the Inventory Control Points 
(ICPs) in managing repairables, and a description of the 
mathematical models in use today for the management of 
repairables at one of the Navy's ICPs, the Ships Parts 
Control Center (SPCC). This is followed by a review of two 
multi-echelon inventory models developed specifically for 
the military: (1) the Multi-Echelon Technique for 

Recoverable Item Control (METRIC) developed by Rand 
Corporation in 1966 for the Air Force; and (2) the 

Availability Centered Inventory Model (ACIM) developed by 
CACI in 1981 for the Navy's use in determining consumer 
level stockage quantities for selected equipments. Both 
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models recognize that the purpose of a supply system is to 
provide sufficient support so that a weapon system is opera- 
tional when it is needed. 

The objective of ACIM is to determine stock levels for 
all repair parts in the equipment in addition to the stock 
levels for the repairable item, while considering where each 
item should be stocked (i.e., what echelon), such that the 
Mean Supply Response Time is minimized subject to a given 
inventory budget. Mean Supply Response Time (MSRT) is the 
mean time it takes the supply system to respond to the 
demand for a replacement part or component. The current 
Navy supply-system goal is 125 hours for ships in CONUS and 
135 hours for ships EXCONUS [Ref. 4: Ch. 4]. ACIM assumes a 
policy of one-for-one ordering between echelons and 
repairing at the wholesale level. Attrition is assumed to 
not occur. Finally, it should be noted that both the ACIM 
and METRIC models are extremely difficult to use because of 
the level of detail of the data and the long computational 
times . 

The Apple Model is a multi- echelon model designed for 
the Navy's wholesale level of repairables and also focuses 
on Mean Supply Response Time, while including a specified 
protection level at the shipboard level as an input param- 
eter. It shows that unless economic reasons dictate, the 
supply system should follow a one-for-one ordering policy 
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for stock lost through attrition and should have a one-for- 
one repair policy also; i.e., no batching of procurements or 
repairs . 

The protection level specified at the shipboard level 
consists of stocks of repairable components which are main- 
tained on board to repair weapon systems while ships are 
deployed and without access to the wholesale supply system. 
By varying this protection level, Apple was able to identify 
situations where the Navy's MSRT goal could not be obtained 
because the shipboard protection level was inadequate. 

The following are the assumptions applicable to his 
model : 

1 failures are generated by a Poisson process; 

2 ships use a one-for-one reorder policy for stock 
authorized on board; 

3 the minimum protection level of spares is the same for 
all ships; 

4 designated overhaul points (DOPs) are established for 
all items; 

5 attrition of items, due to not being turned in or being 
beyond economic repair, is allowed; 

6 repair batch size and procurement lot size are input 
parameters which are determined outside the model; 

7 all demands for stock are satisfied by the wholesale 
system - no lateral resupply (i.e. no resupply between 
ships or NSCs ) ; and 

8 times used throughout are average times expressed in 
quarters . 

Analyses of Apple's model showed that: 
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the wholesale stock level of an item can be greatly 
affected by the repair and procurement policies in 
effect; i.e., a lower stock level is required, if one- 
for-one repair and procurement policies are used. Also, 
by reducing the repair time required for an item, the 
stock level can be reduced ... and by increasing the ship- 
board stock level in our computations, the stock 
required by the wholesale system was greatly reduced.... 
[Ref. 3: p. 86] 



C. TWO RESUPPLY CYCLES 

Two resupply cycles exist in the U. S. Navy - the repair 
and the procurement cycles. Each is discussed separately. 

The repair cycle resupplies the wholesale system by 
receiving not-ready- for issue (NRFI) carcasses from ships, 
inducting them into the repair cycles at a DOP , repairing 
the carcasses, and returning them to the wholesale system in 
a ready-for-issue (RFI) state. The times that affect the 
turnaround time in the repair cycle are: Tl, T2 , RTAT , T3 , 
and any delay resulting from batching of repairs. 

Utilizing Ross [Ref. 5: p. 152], Apple derives the 

average time added to the repair cycle, W(R) , given that the 
repair batch size has been predetermined to be R. His 
formula is presented as equation 2.1 (Note: When batch size 
R equals 1, W(R) =0.) 

W(R) = (R-1)/(2*D*RSR*CRR), (eqn 2.1) 

where: R: repair quantity; 

D: quarterly total expected failure rate; 

RSR: repair survival rate - rate at which NRFI 
carcasses survive the repair process and are 
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returned to RFI condition; 



CRR: carcass return rate - rate at which NRFI 
carcasses are returned to the wholesale 
system from the ship for induction into the 
repair process . 

Combining all times that affect the turnaround time in the 
repair cycle, the mean length of the repair cycle (TT1) is: 

TTl = T1 + T2 + RTAT + T3 + W(R) (eqn 2.2) 

or 

TTl = CRT + RTAT + ( (R- 1) / (2*D*RSR*CRR) ) , (eqn 2.3) 

if we let CRT be the carcass return time (equal to the sum 
of T1 and T2 ) . 

The procurement cycle replenishes the wholesale system 
by procuring new items to replace those which have attrited. 
Attrition occurs in the system due to items not being turned 
in for repair or being beyond economic repair. The times 
that affect the mean procurement time in the procurement 
cycle are: T5 , ALT, PLT , T3 , and any delays resulting from 

the batching of attrited units to accumulate an economic 
order quantity before placing a procurement order. 

Similar to the procedure used above for the repair 
cycle, Apple calculates the delay from accumulating attrited 
units into a batch size of Q before procuring to be: 

W(Q) = ( Q- 1 ) / ( 2*D* ( 1 - ( RSR*CRR ) ) ) . (eqn 2.4) 
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Again, note that when Q equals 1, W(Q) = 0. The mean length 
of the procurement cycle (TT2) is obtained by combining the 
pertinent time variables to obtain: 

TT2 = T5 + ALT + PLT + T3 + W(Q) (eqn 2.5) 

or 

TT2 = PCLT+T5+ ( ( Q- 1 ) / (2*D* (1- (RSR*CRR) ) ) ) (eqn 2.6) 



where: PCLT : procurement cycle lead time (equal to ALT + 

PLT + T3). 

With both the mean repair cycle and the mean procurement 
cycle times known, Apple then develops the mean resupply 
time and the mean number of units in resupply. 

The mean resupply cycle time is obtained by multiplying 
the mean length of the repair cycle (equation 2.3) by the 
probability that a failed unit can be returned to RFI condi- 
tion through repair (RSR*CRR) and summing this with the 
product of the mean length of the procurement cycle (equa- 
tion 2.6) and the probability that an item must be replaced 
through procurement ( 1- (RSR-CRR) ) . Thus, the equation for 
the mean resupply cycle time is: 

MU = ( RSR-CRR ) * ( CRT + RTAT + (eqn 2.7) 

( (R-l)/ (2*D*RSR*CRR) ) ) + ( 1- (RSR*CRR) ) 

*(PCLT + T5 + (Q- 1 ) / (2*D*(1- (RSR*CRR) ) ) ) 
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The mean number of units of item i in resupply i) is 

simply the demand or failure rate (Di) times the mean 
resupply cycle time (equation 2.7) which yields: 

pi = Di * MUi (eqn 2.8) 

D. THE OBJECTIVE FUNCTION 

The expected number of backorders at a randomly selected 
time is equivalent computationally to the total expected 
time-weighted units short (TWUS) per unit of time [Ref. 6: 
p . 185] . That is : 

TWUSi(SWi) = ( pi - SWi ) + (eqn 2.9) 

SWi-1 

E (SWi - xi)*pi(xi; pi). 
xi = 0 

By dividing the TWUSi by the total expected failure rate, 
Di , Apple obtained the average delay per failure or the mean 
supply response time for item 'i'. Adding this to T4, which 
accounts for shipping time from the NSC to a given ship, 
gives the mean supply response time of the wholesale system 
for a given ship. This is expressed as: 

MSRTi = T4 + MSRTRSi ( SWi ) , (eqn 2.10) 

where: MSRTRSi (SWi ) : mean supply response time for the 

resupply cycle (which is equivalent to 
TWUSi (SWi)/Di) . 
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He then derives the average MSRT across all ships for 



item 'i' to be: 



J 

MSRTi(SWi) = E Bi j ( SSi j , SWi ; 0 ij)/Di, (eqn 2.11) 

j=l 



where : Bi j : 

SSij : 

0ij : 



SWi: 



expected number of backorders for item 'i' at 
a random selected time for ship 'j'; 
ship stock level for item 'i'; 

MSRTi-'Dij is the mean demand at ship 'j' for 
item 'i' during an average resupply time and 
is a function of SWi; 

Wholesale stock level of item ’ i ’ ; and 



J 

Di : E Di j . 

j=l 

MSRTi(SWi) is constrained to be no larger than the MSRT 



goal . 

The average supply system MSRT over all items then 
becomes : 



I I 

MSRT = E Di*MSRTi ( SWi ) / E Di . (eqn 2.12) 

i=l i=l 

Through an iterative process, Apple's model can be used to 
add units of stock at the wholesale level until the MSRT 
reaches the specified MSRT goal. 
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The Navy consistently faces funding limitations which 
necessitate prioritization of requirements. This is accom- 
plished by identifying weapon systems with respect to their 
criticality or essentiality. The Navy is able then to deter- 
mine which items will reduce the difference between what is 
required and what is available given a budget while doing 
the least amount of damage to the established MSRT goal. 
Apple incorporated this concept of essentiality in his model 
and then stated the problem as: 

I 

Minimize: 2 Ei*Di*MSRTi ( SWi ) 

i=l 

I 

Subject to: 2 Ci "'SWi < B, 

i = l 

where: Ei : Item Mission Essentiality Code associated with 

item "i"; [Ref. 4: p. 4-40] 

Ci: procurement cost of item ’ i ’ ; 

B: total budget dollar constraint. 

Although Apple's problem was presented as above, he 
actually had initially stated it as: 

The objective is to find the level of wholesale stock, 
SWi, (consisting of both RFI and NRFI assets) for each 
of the 'i' items in the supply system required either to 
minimize the MSRT subject to a budget constraint or to 
determine the minimum cost solution which attains a 
predetermined MSRT goal. [Ref. 3: p. 60] 
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Apple analyzed the former. The latter can be written as 
follows : 



I 

Minimize: £ Ci*SWi 

i=l 

I I 

Subject to: £ Ei*Di*MSRTi ( SWi ) / £ Ei*Di < MSRT goal 

i=l i=l 

This problem statement will be the basis for the aggregate 
demand model which will be developed in Chapter III and 
analyzed in Chapter IV. 

The final step of Apple's thesis was to use marginal 
analysis in the selection of a wholesale stock level to meet 
a specified MSRT goal at the shipboard level for an example 
system consisting of one ship (or customer). 
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III. THE AGGREGATE DEMAND MODEL 



This chapter will develop a conversion of Apple's model 
to an aggregate demand model which retains the same assump- 
tions as outlined in Chapter II. 

A. AGGREGATION OF DEMAND 

The Navy's ICPs are concerned with managing the whole- 
sale supply system based on an aggregate demand from many 
customers (i.e., many ships, shore stations. Foreign 
Military Sales, other branches of the Department of Defense, 
etc.). This aggregate demand is currently forecasted in the 
Uniform Inventory Control Programs (UICP) , which are various 
computer programs that the Fleet Material Support Office 
(FMSO) has developed to provide the ICPs with scientific 
inventory management techniques. Since this capability is 
available, it is appropriate to utilize the forecasted 
aggregate demand as an estimate of Di . By using this aggre- 
gated data, the problems of obtaining customer level data, 
such as Apple's model requires, are avoided. 

B. OBJECTIVE FUNCTION AND MSRT GOAL 

If we follow the current UICP approach, we would set an 
MSRT goal for the supply system just as the current Supply 
Material Availability goal of 85% has been set by the Chief 
of Naval Operations. Selection of that goal implies a 
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wholesale stock (SWi) level for each item (i) and its 
related monetary value or investment level in the system. 
Summing these investment levels, we get: 

I 

L Cli*SWi 
i = l 

where: Cli: unit procurement cost or price. 

Ideally, it would be nice to always possess the 
resources necessary to achieve the lowest possible MSRT 
(i.e. zero, where all items are always immediately avail- 
able). However, in reality, the Navy has funding limitations 
that restrict the size of the Navy Stock Fund (NSF). The NSF 
is a revolving fund managed by the Naval Supply Systems 
Command and consists of money and/or stock. The stock fund 
is reimbursed by the customer when stock is issued and these 
funds are used to procure new items or to repair NRFI items 
to replace the inventory that has been issued. Therefore, a 
continual concern for the Navy is that investment in whole- 
sale stock be minimized to achieve the established MSRT 
goal . 

Because of our concern over the money tied up in invest- 
ments, we would like to minimize it. Therefore, we can write 
our problem now as: 

Find SWi which 
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I 

Minimizes: E Cli*SWi 

i = l 

I I 

Subject to: E Di*Ei*MSRTi ( SWi ) / E Di*Ei << MSRT Goal 

i=l i= 1 

where: MSRTi(SWi): the mean supply response time for the 

wholesale system for item (i) when the 
wholesale stock level of item (i) is 
SWi. 

C. ITEM MISSION ESSENTIALITY CODE 

Recent analyses have developed a replacement for Ei 
known as the Item Mission Essentiality Code (IMEC). Each 
item (i) is categorized by this code which is a weighting 
factor to indicate that certain weapon systems are more 
critical than others. IMECs are defined in NAVSUP 
Publication 553 [Ref. 4: p. 4-40] as follows: 

IMEC Definition 

4 Loss of primary mission capability 

3 Severe degradation of a primary mission 
capability 

2 Loss of a secondary mission capability 
1 • Minor mission impact 

They were chosen to replace the Ei values for Apple's model. 
The problem with using these integer values, as Apple 
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observed, is Chat an item with an IMEC of 4 is not neces- 
sarily twice as important as an item with an IMEC of 2. 
This problem can be overcome by establishing four separate 
MSRT goals to correspond to the four IMEC levels and by the 
deletion of the IMEC in the MSRT constraint. 

The aggregate model will thus assume the separation of 
all items into IMEC categories and the assignment of an 
appropriate MSRT goal for each level (i.e. possibly 24 hours 
for IMEC 4; 72 hours for IMEC 3; 125 hours for IMEC 2; and 
150 hours for IMEC 1). 

D. THE AGGREGATE MODEL SOLUTION 

Our problem now is to find SWi for all i = 1,2,..., I, 
(having a given IMEC code) which 



I 

Minimizes: I Cli*SWi 

i=l 



I I 

Subject to: E Di'-MSRTi (SWi) / E Di < MSRT Goal. 

i=l i= 1 



In order to compute the aggregate MSRTi(SWi) value 
corresponding to a given level of wholesale system stock, 
SWi, it is necessary to recall equation 2.10 from Chapter 
II:- 

MSRTi(SWi) = T4 + MSRTRSi (SWi ) (eqn 3.1) 
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For simplification we will assume T4 to be zero since 
shipment of a RFI unit to a customer can be expected to take 
negligible time relative to RTAT and PCLT. Thus, the aggre- 
gate MSRTi may be rewritten as: 

MSRTi(SWi) = TWUSi ( SWi ) /Di (eqn 3.2) 

and the constraint is then 
I I 

Z TWUSi (SWi)/ Z Di < MSRT Goal . (eqn 3.3) 

i=l i=l 

An iterative process can be utilized to search for each 
optimum SWi. The initial step in the solution of this 
problem is to obtain the mean number of units in resupply 
(jLti) which was derived in Chapter II and resulted in equa- 
tion 2.8. Expanded, this equation becomes equation 3.4 where 
the argument i has been surpressed for clarity. 

V = D* ( ( RSR “ CRR ) * ( CRT + RTAT + (eqn 3.4) 

((R - 1)/ (2*D*RSR*CRR) ) ) + ( 1- (RSR--CRR) ) 

*(PCLT+T5 + ( Q - 1 ) / ( 2'- v D“ ( 1 - ( RSR*CRR ) ) ) ) ) 

Values for all variables in equation 3.4 are available 
from forecasted and historical data maintained in the UICP 
files except for repair induction (R) and procurement (Q) 
quantities. These two variables are specified as a result of 
budgetary and policy considerations. For this model, both Q 
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and R will equal 1 in the base model which is presented in 
Chapter IV and they will be allowed to vary in the subseq- 
uent parametric analyses. 

Once y i is known, it is easy to calculate the total 
expected time-weighted units short (TWUS) per quarter for 
each item (presented in equation 2.9 and repeated here for 
convenience ) . 

TWUSi(SWi) = ( pi - SWi ) + (eqn 3.5) 
SWi-1 

E (SWi - xi)*Pi(xi; pi). 
xi=0 

We next calculate TWUSi for SWi = 0 for all items. We then 
combine the result with each item's forecasted demand, as 
shown in the left-hand side of inequality 3.3, to arrive at 
the system-wide MSRT that will be provided when SWi equals 
zero for all items. This calculated MSRT, denoted as CMSRT , 
is compared to our MSRT goal and, if CMSRT is less than or 
equal to MSRT, we stop with SWi being zero across all items. 

If CMSRT when all SWi = 0 is greater than the MSRT goal, 
we implement a marginal analysis procedure to determine SWi. 
This procedure makes use of a weighting factor for each item 
of stock based upon cost and time-weighted units short. 
This expression is represented as: 

WTi = Cli/ (TWUS (SWi- 1 ) - TWUS(SWi)). (eqn 3.6) 
This ratio expresses the increase in investment cost of each 
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item relative to the benefit in reduced response time 
derived from adding one additional unit of the item to the 
wholesale stock. 

For each item being considered, we compute WTi assuming 

SWi = 1 and then add one unit to that item k for which WTk = 

min { WTi } . We again check to see if the MSRT goal is 
i 

satisfied by computing the left-hand side of the constraint 
(inequality 3.3) and comparing it to the MSRT goal value. If 
the computed MSRT is greater than the MSRT goal, a new value 
of WTk is computed assuming SWk = 2 before comparing it with 
other WT values. Again, we select that item having the 
smallest WTi and increase its wholesale level by one unit. 
This process continues until the computed MSRT is less than 
or equal to the MSRT goal. 

Finally, with all SWi values known from this last step 
of the marginal analysis procedure, the value of the objec- 
tive function is computed by summing Cli - SWi over all items. 
This will provide approximately the minimum total investment 
required to meet the given MSRT goal. 

The SWi values calculated by the model represent the 
maximum values of the inventory position. As demands occur 
the inventory position will decrease. When a repair induc- 
tion is made, the inventory position for item i will be 
increased by the value of the expected successful regenera- 
tions (or Ri/RSRi). When the inventory position immediately 
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after an induction reaches or falls below SWi - Qi , a 
procurement of Qi should be made. This will immediately 
return the inventory position to SWi. 

E. PROCUREMENT AND REPAIR INDUCTION VALUES 

The two variables playing significant roles in both this 
aggregate model and Apple's model are the procurement quan- 
tity Qi and the repair induction quantity Ri . Their values 
have a major impact upon the wholesale stock levels (SWi) 
required to achieve the specified MSRT goal. Therefore, a 
brief discussion concerning Qi and Ri and their relationship 
to the model and to existing budgetary practice in the Navy 
is appropriate prior to discussing the model analysis and 
results . 

Apple considered a problem for making a one-time buy of 
SWi, given quantities of size Qi and Ri would be bought and 
inducted for repairs, respectively, whenever necessary. In 
reality, this may not be possible. The annual UICP budgets 
are designed to pay for procurements and repairs (rather 
than buying SWi), but the amounts received from Congress may 
not be sufficient to buy Qi and Ri whenever necessary. Thus 
a limited UICP budget imposes additional constraints on our 
problem. These are: 

I 

Z (Cli * ni * Qi) BP (procurement budget) ; (eqn 3.7) 
i = l 
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I 

£ (C2i * mi * Ri ) < BR (regeneration budget) (eqn 3.8) 

i = l 



Cl: 


unit procurement 


cost 


or price; 


C2 : 


unit repair cost; 






n : 


number of procurement 


buys per year; 


m: 


number of repair 


inductions per year; 



The annual number of buys and the annual number of 
inductions would depend on the expected annual number of 
attritions and regenerations. Thus: 

n = Number of attritions/Q (eqn 3.9) 

and 

m = Number of regenerations/R . (eqn 3.10) 

When we introduce these expressions into the budget 
constraints above, we get: 

I 

£ (Cli * Number of i attritions)^ BP (eqn 3.11) 

i=l 



I 

£ (C2i * Number of i regenerations) < BR (eqn 3.12) 

i=l 



An implicit assumption of the Apple model is that the 
annual number of attritions can always be bought and the 
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annual number of regenerations can always be funded. If BP 
and BR do not allow this, then the wholesale stocks cannot 
sustain the SWi levels and the MSRT provided by those 
levels . 

The problem of a limited budget is currently handled by 
adjusting the reorder points of the current UICP models each 
year. In particular, the reorder points are lowered and 
result in postponement of buys and inductions when the 
budget constraints are severe. Thus, Qi and Ri are not 
necessarily always procured or inducted for repair "whenever 
necessary". Any new repairables model must therefore have 
provisions for handling these severe budget constraints. 

The incorporation of the constraints, given by inequali- 
ties (3.11) and (3.12), into the aggregate model would make 
the model much more complex. Thus, this thesis will not 
attempt to solve this larger problem. Instead it will assume 
Qi and Ri to be given, derived perhaps from some initial 
optimization step. Once parametric analyses have been 
conducted and the impact of Qi and Ri are well understood, a 
methodology for optimization of the three constraint problem 
may become obvious. 
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IV. MODEL RESULTS AND PARAMETRIC ANALYSIS 



This chapter begins with a brief description of the 
computer program developed for solving the aggregate demand 
model discussed in Chapter III. This description includes 
the characteristics of the input and output variables. The 
rest of the chapter is devoted to several parametric anal- 
yses . 

A . COMPUTER PROGRAM 

A flow chart of the solution procedure for the aggregate 
demand repairables model derived in Chapter III is provided 
in Appendix A. The procedure was programmed in FORTRAN and 
run with the WATFIV compiler on the IBM 3033 at the Naval 
Postgraduate School. Appendix B contains a listing of the 
program . 

The program's input and output variables are detailed in 
Table I. As was mentioned in Chapter III, MSRT represents 
the mean supply response time goal of the wholesale system 
for a group of items having the same Item Material 
Essentiality Code (IMEC). As a consequence, there is no need 
to incorporate an essentiality weighting factor in any of 
the formulas used in this program. CMSRT is the computed 
mean supply response time for a given set of input variable 
values. Both CMSRT and MSRT are stated in days in the 
output . 
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TABLE I 



VARIABLE DESCRIPTIONS 



Input Variables 
N: Number of line items 

MSRT : Mean Supply Response Time Goal (in quarters) 

NUN : National Item Identification Number 

ATT: Carcass Attrition Rate (in number of units 

per quarter) 

REG: Carcass Regeneration Rate (in number of units 

per quarter) 

RSR: Carcass Repair Survival Rate (probability) 

CRR: Carcass Return Rate (probability) 

CRT: Carcass Return Time 

RTAT : Repair Turn Around Time 

PCLT : Procurement Cycle Lead Time 

D: Quarterly demand rate 

Q: Procurement size 

R: Repair induction size 

Cl: Cost to procure one unit 

C2: Cost to repair one unit 



Output Variables 

CMSRT : Computed mean supply response time (in days) 

MSRT: Mean supply response time (in days) 

SW: Computed wholesale stock level 

COSTSW: Total investment cost per line item 

INVEST: Minimum total investment required over all 

line items 



As with Apple's model, the probability distribution 
assumed for demand is the Poisson. For those items with a 
mean number of units in resupply (fii) greater than 20, a 
normal distribution with continuity corrections was used as 
an approximation. 
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B. PARAMETRIC ANALYSES 



1 . Base Case Data Set 

The parametric analyses of the proposed aggregate 
demand model used a reference or base set of data to facili- 
tate analytic comparison when input parameters were varied. 
The complexity of this model, with respect to the number of 
parameters, both for the system (i.e., R, Q, and MSRT 
values) and for each individual item (i.e., RTAT, CRR , CRT, 
RSR, PCLT, Cl, and D) , required that the base case also fix 
as many factors as feasible at common values across all 
items. In fact, all were so fixed except for the unit costs 
Cl and C2. In addition, the number of items was limited to 
two. The common input data values for the two items are 
provided in Appendix C. 

The key system input variables in the model are the 
procurement and repair induction quantities, Q and R respec- 
tively, and the Mean Supply Response Time (MSRT) goal. 
Although Apple [Ref. 3] argued that the optimum values for Q 
and R are unity, larger values of Q and R are often neces- 
sary because of budgetary constraints. However, as a base 
case for reference, both Q and R are assumed to be unity. 

The MSRT goal for the base case was 125 hours (or 
.0572 quarters) which is the same goal established for CONUS 
ships [Ref. 4], and discussed in Chapter II. 

Thus, the base case consisted of a two-item popula- 
tion, with neither item batched for repair or procurement 
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R=Q=1). The only varying input parameters were the 



(i.e. , 

unit costs Cl and C2. For item one Cl was $10,000 and C2 
was $5,000 and for item two Cl was $150,000 and C2 was 
$75,000. 

For the base case, the minimum investment levels of 
wholesale stock were 33 and 27 units for items one and two, 
respectively. The total aggregate dollar investment was 4.38 
million dollars. The computed system MSRT was 4.7224 days 
which was the closest the solution could come to the goal 
(5.2052 days) without exceeding it. 

2. Effects of Varyin g Repair Induction Quantities 

As discussed in the previous section the base case 
assumed no batching for repair or procurement. However, 
because a severe budget constraint may require batching, it 
is important to study the impact of batching. Although 
batching of both Q and R quantities is feasible, our anal- 
ysis considers only the results of variations to the repair 
induction quantity, R. For our sample input data the model 
was insensitive to large changes in Q. 

Table II provides, in matrix form, the results of 
thirty-six combinations of R for the two items, denoted by 
NUN 1 and NUN 2. For ease of discussion, each cell in the 
matrix is identified by (row, column) . Contained within each 
cell are the quantities of wholesale stock (SW) calculated 
from the model assuming the specified R1 and R2 values and 
the established MSRT goal of 5.2052 days. The SW1 for NUN 1 
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is in the upper right of each cell and SW2 for NUN 2 is in 
the lower left. The base case results are found in cell 
( 1 , 1 ). 

The impact on investment level when the R2 quanti- 
ties for NUN 2 are held constant and the R1 values for NIIN 
1 are allowed to vary from one to six can be seen by consid- 
ering each row. The SW2 level for NIIN 2 remains constant 
for a specified R2 value, while the SW1 level for NIIN 1 
increases from either 33 or 34 to a high of 36, as the batch 
size, Rl of NIIN 1, increases. This analysis suggests that 
in the aggregate demand model, as the R value for a single 
item increases, the investment level for that item either 
remains constant between cells or it increases by one unit. 
This is graphically depicted in Figure 4.1 and Figure 4.2. 
Note that Figure 4.1 corresponds to the odd values of R2 , 
and Figure 4.2 corresponds to the even values. The differ- 
ences between the two figures are addressed below. 

The R2 analysis obtained similar results as the Rl 
analysis. The results of holding Rl values constant and 
allowing the R2 values to vary from one to six can be seen 
by considering each column. Figure 4.3 shows the stepwise 
behavior (it is the same for all Rl values). 

An interesting interaction effect can be observed in 
column one, where the SW1 value alternated between 33 and 34 
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TABLE II 



INVESTMENT LEVELS (SW1.SW2) 
AS A FUNCTION OF R1 AND R2 



R2 



R1 



33 \ 34 \ 35 \ 35 \ 36 \ 36 

27 \ 27 \ 27 \ 27 \ 27\ 27 



34 \ 34 \ 35 \ 35 \ 36 \ 36 

28 \ 28 \ 28 \ 28 \ 28\ 28 



33 \ 34 \ 35 \ 35 \ 36 \ 36 

28 \ 28 \ 28 \ 28 \ 28\ 28 



34 \ 34 \ 35 \ 35 \ 36 \ 36 

29 \ 29 \ 29 \ 29\ 29\ 29 



33 \ 34 \ 35 \ 35 \ 36 \ 36 

29 \ 29 \ 29 \ 29 \ 29\ 29\ 



34 \ 34 \ 35 \ 35 \ 36 \ 36 

30 \ 30 \ 30\ 30\ 30\ 30 
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even though the Rl remained constant at one. The reason for 
this oscillation is not obvious and further study of this 
situation is needed. 

Several other results are worthy of comment. First, 
Table III is a matrix presentation of the computed MSRT (in 
days) where the model terminated once the MSRT goal of 
5.2052 days was achieved. These computed MSRT figures 
consistently remained constant across rows. This is attrib- 
uted to the fact that the model stops increasing the invest- 
ment level of SW1 when the computed time weighted units 
short (CTWUS) is either zero or so close to zero that NUN 1 
does not contribute any significant value to the numerator 
of equation 3.3 (used to calculate MSRT). The model then 
continues to add to the investment level of SW2 until the 
MSRT goal is achieved. Actually, it turned out that the 
normal approximation gave small negative values for CTWUS1. 
These continued to be used for the marginal analysis but the 
MSRT values were assumed to be zero as soon as CTWUSl went 
negative. The reason for the normal approximation giving a 
negative CTWUSl is not clear and needs to be studied 
further. Perhaps the break point for using the normal 
approximation should be much larger than the value of 20 
assumed for the aggregate demand rate. 

The column behavior of CMSRT is different from the 
row behavior. This is because the CTWUS for NUN 2 does not 
reach a value near zero until R2 is much larger than six. 
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Figure 4.1 SW1 Sensitivity to R1 when R2 = 1. 
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Figure 4.2 SW1 Sensitivity to R1 when R2 = 2. 
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Figure 4.3 SW2 Sensitivity to R2 for 1 >< R1 < 6. 
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COMPUTED MSRT 


TABLE III 

WITH REPAIR INDUCTION QUANTITY 


1 

VARIED 




1 


2 


R1 

3 


4 


5 


6 


1 


4.7224 


4.7224 


4.7224 


4.7224 


4.7224 


4.7224 


2 


3 . 9546 


3.9546 


3.9546 


3.9546 


3 . 9546 


3.9546 


3 


4.7980 


4.7980 


4.7980 


4.7980 


4.7980 


4.7980 


R2 

4 


4.0280 


4.0280 


4.0280 


4.0280 


4.0280 


4.0280 


5 


4.8559 


4.8559 


4.8559 


4.8559 


4.8559 


4.8559 


6 


4.0983 


4.0983 


4.0983 


4.0983 


4.0983 


4.0983 



Comparison of the matrices of Table II and Table III 
shows that two consecutive cells, in a column may have the 
same SW for their respective NIINs, but the computed MSRT 
(CMSRT) varies. An example of this occurs in cells (2,3) and 
(3,3), where SW1 and SW2 remain at 35 and 28, respectively, 
while the CMSRT in cell (2,3) is 3.9546 days and 4.7908 days 
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in cell (3,3). The explanation for this rests with the 
manner in which the value of R2 contributes to the calcula- 
tion of time weighted units short (CTWUS). As a result of R2 
increasing by one unit the CTWUS for NUN 2 is increased 
slightly, which, in turn, causes the CMSRT to increase. 
However, the increase in CMSRT is not large enough to exceed 
the MSRT goal . If the CMSRT had exceeded the MSRT goal , an 
additional unit of NUN 2 would have been added to SW2 (This 
is what happens in cell (4,3)). The reason that no similar 
change is observed when we examine neighboring row cells is 
that CTWUS1 is essentially zero as discussed earlier. 

The final observation with respect to these matrices 
is that the total investment level SW1 + SW2 along the main 
diagonal always increased between cells. This continues 
until both SW levels reach a point where any further addi- 
tions would not provide any improvement to the calculated 
MSRT since each CTWUS would be zero. In this example, as 
discussed earlier, this level would be where SW1 and SW2 
equaled 33 and 36 respectively. 

3 . Wholesale Investment Stock Level Impact on MSRT 

Because the Mean Supply Response Time (MSRT) goal is 
one of the major system input parameters, we should also 
study the impact of changing this parameter. To accomplish 
this efficiently, two cells were selected arbitrarily from 
Table II (cells (1,6) and (6,1)), and the MSRT goal was set 
at zero days. Then one wholesale stock level (SW) was 
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allowed to increase while the second SW was held fixed. 
With each addition to SW, a new CMSRT was computed and 
represents the attainable MSRT for that step. (The basic 
computer program for the model, presented in Appendix B, was 
modified to conduct this analysis. The modifications are 
provided in Appendix D). 

Figure 4.4 shows the results for cell (1,6). This 
shows how CMSRT decreases as SWl increases when R1 = 6 , R2 = 
1, and SW2 =27. A similar result is obtained (Figure 4.5) 
when Rl = lj R2 = 6, SWl = 34, and SW2 is allowed to 
increase . 

This analysis confirms an observation from the 
previous section that beyond a certain SW value for each 
item the computed time weighted units short (CTWUS) is 
essentially zero and additional investment in that SW will 
not result in an improvement in CMSRT. As a consequence, it 
seems appropriate to set a bound on each item's CMSRT so 
that further investments won't be made. Reference 2 found a 
bound of 0.001 days to be reasonable. 

C. CONCLUSIONS 

This chapter described the computer program used to 
compute optimal SWi and has presented preliminary parametric 
analyses of the aggregate demand repairables model proposed 
in Chapter III. Based upon this limited analysis, some 
conclusions may be drawn. 



47 



o 




WHOLESALE STOCK LEVEL 



Figure 4.4 Sensitivity of Attained MSRT to Changes SW1 
when Rl =6, R2 = 1, and SW2 = 27. 
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Although Apple showed that optimal Q and R are equal to 
one, we realize that the values of Q and R may not be able 
to be that small and may be larger than unity for some items 
because of severe budget constraints. Their values will 
depend on the relative unit costs of the various items. 

The model also assumes that Q and R are bought and 
inducted, respectively, whenever necessary. This implies 
that annual procurement and repair budgets must be adequate 
to fund these procurements and repairs once Q and R are 
selected. If later years' budgets are less than that used to 
determine Q and R, then quantities of sizes Q and R cannot 
continue to be bought and repaired. Larger Q and R values 
are then needed so that the total annual buys and repairs 
cost less. The result is that a larger investment in whole- 
sale stock levels and hence an increase in the stock fund 
ceiling will be required to achieve the established MSRT 
goal or the goal must be lowered. 



49 



o 

m 




NHOLESRLE STOCK LEVEL 



Figure 4.5 Sensitivity of Attained MSRT to Changes in SW2 
when R1 = 1, R2 = 6, and SW2 = 34. 
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V. SUMMARY AND RECOMMENDATIONS 



A . SUMMARY 

A brief overview of the repairables system was provided 
as an introduction to the review of Apple's thesis [Ref. 3] 
in Chapter II. His model is a performance and Navy oriented 
multi-echelon model using mean supply response time (MSRT) 
as an objective function. It incorporated considerations of 
both repair and procurement to sustain inventories at the 
wholesale and shipboard levels. It was identified as being 
less complex from a computational perspective than other 
performance oriented models (i.e., ACIM and METRIC). 
Nonetheless, it was clear from its development that many 
parameters must interact to obtain a solution to the problem 
of minimizing MSRT. 

Chapter III used Apple's model as a foundation to 
develop an aggregate demand model which requires fewer 
parameters than Apple's and can use existing data from the 
UICP. The objective of this model was to minimize the 
wholesale stock level investment while attempting to achieve 
an established MSRT goal. By stating the model in this way, 
it is able to address both the issue of readiness, which is 
implied by the MSRT goal, and the Navy's concern over 
investment levels since they are directly related to the 
value of the Navy Stock Fund. The Navy must be continually 
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aware of its position with respect to the Navy Stock Fund so 
that stock levels are not built up unnecessarily when it 
does not aid in achieving the desired MSRT goal. 

Finally, grouping of items with the same IMEC and estab- 
lishing individual MSRT goals for each IMEC category allowed 
the model to account for essentiality. Marginal analysis was 
proposed as the optimization technique for this model. 

The last issues addressed in Chapter III were the annual 
procurement and repair budget constraints. It was noted that 
when these two constraints were added then marginal analysis 
could no longer be used for optimization. No alternative 
optimization technique has been developed as yet to solve 
this larger problem. 

Chapter IV presented limited parametric analyses of the 
aggregate model for the case of two items and assumed iden- 
tical values for most of the various item parameters. The 
procurement costs, repair costs, and repair quantities were 
allowed to differ. In particulular , the effect of varying 
the repair induction quantity was examined. The results 
showed that the required wholesale stock investment level 
had to be increased when the repair quantity is increased in 
order to achieve an established MSRT goal. 

Chapter IV also showed that by holding all parameters 
constant and allowing one item's SW level to increase, the 
aggregate attained MSRT decreased exponentially. A point of 
diminishing return is reached such that continuing to 
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increase an item's investment level results in inefficient 
use of resources. A lower bound on a given item's MSRT was 
proposed to prevent such inefficiency. 

B. RECOMMENDATIONS 

The most critical area requiring further research is how 
to incorporate the annual procurement and budget constraints 
into the model. Annual procurement and regeneration budgets 
typically change yearly. Implementation of the aggregate 
demand repairables model by an ICP requires that these 
varying budget constraints be accommodated. In addition, the 
impact of these constraints on attainable MSRT goals and 
investment levels must be understood. 

A procedure must be developed for determining the 
optimum procurement (Q) and repair induction (R) sizes which 
will be feasible given the annual procurement and regenera- 
tion budgets. If the annual budgets are sufficient to fund 
all procurements and regenerations when Q and R are unity, 
then obviously both should remain at unity for input to the 
model since the investment levels will be lowest. If the 
budgets are insufficient to allow this, larger values of Q 
and R must be determined so that buys and repair inductions 
can be postponed until the next year. Because the maximum 
value of the inventory position must then be increased for 
some items, it is important to increase those which have the 
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best trade-off between the costs of increased investment 



levels and the amounts of procurement and regeneration 
budgets consumed. 

As discussed in Chapter IV, further study is needed to 
determine the reason that the normal approximation gives 
small negative values for CTWUS1. Also, additional study is 
required to identify the cause of the oscillation that 
occurred in the SW1 values of column one in Table II. 

Since this thesis only studied a simplified example, 
further study is needed with a large number of items and 
with changes to the parameters held fixed in Chapter IV. 
Finally, actual data for repairable items from the UICP data 
base should be used in a performance evaluation of the 
aggregate model as compared to the current UICP repairables 
model . 
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APPENDIX A 



AGGREGATE DEMAND MODEL FLOW CHART 
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MSRT : MEAN SUPPLY RESPONSE TIME GOAL 

MU : TOTAL PROGRAM PROBLEM VARIABLE 

MUP : PROCUREMENT PROBLEM VARIABLE 

MUR : REPAIR PROBLEM VARIABLE 

NUN : NATIONAL ITEM IDENTIFICATION NUMBER 
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REG ( I ) =ARR (1,2) 
Cl (I) =ARR (1,3) 
C2 ( I ) =ARR (1,4) 
RSR ( I ) =ARR (1,5) 
CRR ( I ) =ARR (1,6) 
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*CMSRT,Q(N) ,R(N) 

CALL CTWSWO (N , SW , CTWUS1 , MU , P , C ,CMSRT ,MSRT , SWR , D , PCLT ) 
IF (CMSRT.LE.MSRT) GO TO 7090 
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* C1,C2, RSR , CRR , CRT , RTAT , PCLT , D , MU , Q , R ) 

IF (MM .NE.NIIN(I)) GO TO 7075 
GO TO 7085 
7075 J=SWR ( I ) +1 

GO TO 7085 
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7400 SUBROUTINE CPMSRT (CMSRT , N , D , CTWUS 1 ,MSRT) 
INTEGER N 

REAL CMSRT, D(N) ,CTWUS1(N) , MSRT , CCTWUS , DD 
CCTWUS=0 .0 



7410 1 = 1, N 
X=CTWUS1 ( I ) 

IF (CTWUS1 ( I ) . LT . 0 ) 
CCTWUS=CCTWUS + X 
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7600 SUBROUTINE SORTS ( ARR ,N ,NI IN , ATT , REG , Cl , C2 , RSR , CRR ,CRT , 

*RTAT , PCLT , D , WT , SWR , MU , CTWUS 1 , CTWUS 2 , Q , R ) 

INTEGER N,NIIN(N) ,S1,S4, SWR (N) 

REAL ARR (N , 10 ) , ATT(N) ,REG(N) ,C1(N) ,C2(N) ,RSR(N) ,CRR(N) ,CRT(N) , 
*RTAT(N) , PCLT (N) ,D(N) ,WT(N) ,MU(N) ,CTWUS1(N) ,CTWUS2(N) ,Q(N) ,R(N) 



DO 7610 K=JJ,N 

IF (WT(L) .LT.WT(K) ) GO TO 7610 
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NN : INDEX VARIABLE 

MSRT : MEAN SUPPLY RESPONSE TIME GOAL 

MU : TOTAL PROGRAM PROBLEM VARIABLE 

MUP : PROCUREMENT PROBLEM VARIABLE 

MUR : REPAIR PROBLEM VARIABLE 
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C*** ROUTINE TO REARRANGE ARRAY ASSIGNMENTS ******* 
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MUR ( I ) = (RSR ( I ) *CRR ( I ) ) * ( CRT ( I ) +RTAT (I)+((R(I)-1.)/ 
( 2 . *D ( I ) *RSR ( I ) *CRR ( I ) ) ) ) 

MUP(I)=(1.-(RSR(I)*CRR(I) ) )*(PCLT(I)+( (Q(I)-l. ) )/ 
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MU (I)=(D(I)*( MUR ( I ) +MUP ( I ) ) ) 
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